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^ ! We give an overview of the current theoretical status of proper time renormahzation group 

0^ flow equations applied to QCD. These equations give the evolution of coupling constants in 

^ an effective QCD Lagrangian as a function of an infra-red cut-off parameter. This parameter 

Q characterizes the resolution with which the system is looked at. Decreasing resolution trans- 

^ forms partonic quarks into constituent quarks and generates pion bound states. The evolution 

Ph. equations can also be applied to finite temperature and finite density QCD. 

Oh! 

> 

X 



* Supported by the European Network on Electron Scattering off Confined Partons (ESOP) under 
contract no. HPRN-CT-2000-00130 



1 



Quarks and gluons are the ultimate building blocks of hadrons. Their interaction is governed 
by the QCD coupling which grows dynamically for low energies. With the increasing coupling 
gluons and quark-antiquark pairs multiply indefinitely. Therefore QCD forms new structures 
of dominant collective units in the infra-red. The renormalization group is the principal tool to 
investigate the formation of these new structures. Its importance in quantum field theory can 
only be compared to the Schrodinger equation in quantum mechanics. In this paper I use the 
term renormalization group in a wide sense including the evolution of systems of Lagrangians. 
It may well be that the degrees of freedom change under evolution. Such an approach is not so 
well researched yet, but the necessity of a theoretical framework for it - especially in QCD - is 
universally recognized. I will try to give a summary of our work using renormalization group 
techniques to handle quark degrees of freedom evolving into hadronic bound states. 

The quark-hadron transition governs strong interaction physics when the resolution is lowered 
below 1 GeV. Hadronic physics is well researched up to a mass scale of 1.5 GeV. Its spectroscopy 
still lacks a clear identification of glueballs and hybrid states. Yukawa models like the sigma 
model are natural to describe hadron physics at low energies. In order to bridge the gap between 
quark and meson physics we choose a hybrid approach where the original action contains a 
four-quark interaction which is bosonized into sigma and pion fields. These mesonic degrees 
of freedom are not propagating fields in the ultraviolet, i.e., their kinetic term is zero. They 
couple to the quark fields and mimic the local four-fermion interaction. In principle all the 
mesonic couplings arising from the Fierz transformation of a local color exchange interaction 
should be taken into account. Mean field theory shows a basic ambiguity to the possible Fierz 
transformations [|I| . The RG-equations can overcome this problem. At the ultraviolet scale the 
meson-quark Yukawa coupling can be set to unity. There is only the mass parameter of the 
combined sigma and pion fields in the Wigner mode which sets the mass scale for the evolving 
system. The renormalized Yukawa coupling decreases strongly in the infra-red which makes the 
evolution infra-red stable. In the approximation with meson loops we have indications that the 
infra-red evolution has fixed point character. Reversely the divergence of the Yukawa coupling 
for large scales signals the unsatisfactory high energy behaviour of the hadronic theory. Above 
the 1 GeV scale the compositeness of hadronic objects becomes important and one has to choose 
a quark-gluon basis. 

In Fig. 1 we summarize the traditional approach to hadronic and nuclear physics indicating 
the hierarchy between quark physics, nucleon-mesonic physics and nuclear physics. A division 
into these subfields is highly efficient when special topics in a field are researched, because then 
refined methods can be applied to get maximum insight into a detailed aspect. On the other 
hand it is refreshing to cross the borderlines of these fields and see how things are connected. 
Indeed modern high energy experiments in nuclear physics look at nuclei with high resolution 
in electron nucleus collisions and at high excitation energies in relativistic heavy ion collisions. 
They naturally relate nuclear physics to quark physics. 

The renormalization group connects regions of different length scales and is the appropri- 
ate tool for such an interdisciplinary approach. It avoids the appearance of arbitrary cut-off 
functions when additional quantum corrections appear. Going beyond the large Nc appro- 
ximation in the Nambu Jona Lasinio (NJL) model one needs to include meson loops which 
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Figure 1. The traditional separation of quark physics, the physics of the nucleon/delta and 
nucleon-nucleon interaction and nuclear physics 



necessitates an additional cut-off in the Schwinger- Dyson approach. The unifying aspect of 
the renormalization group is shown in Fig. 2 and illustrated by the work here. Starting from 
the renormalization group in the vacuum, which only tests physics along the resolution axis, 
I will cover high temperature physics and high density investigations in the same theoretical 
framework. 
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Figure 2. The renormalization group unifies physics along the three different axis: physics at 
varying resolution, finite temperature and density 



Renormalization group flow equations describe the average of an effective action and represent 
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the continuum analogue of a block spin transformation. I have tried before to use the Monte 
Carlo renormalization group to extend the lattice gluon action of Wilson into the long distance 
region 0,0,0 . This approach introduces new collective variables into the lattice action, namely 
the colordielectric fields. The resulting QCD action has a strong coupling which is modulated by 
the colordielectric field and confines. It reproduces results obtained on the finer lattice with an 
accuracy of 10 %. A large numerical effort is needed to solve a set of overdetermined Schwinger- 
Dyson equations. Already at that time it was desirable to have an analytic scheme to follow the 
flow of couplings when high momenta larger than a coarse-graining scale k are integrated out. 
The flow equations proposed by Wegner and Houghton [§, Polchinsky [0, Wetterich and Berges 
present such a scheme. They are ultraviolet and infra-red finite through the introduction 
of a scale dependent cut-off function. We use renormalization group flow equations with a heat 
kernel cut-off |[T0| , pJ] , p!2[| which we think are more practical. Moreover, they are so simple that 
they are ideally suited to learn and teach renormalization in field theory. No infinities appear at 
any time. All integrals can be done analytically and have a simple form. The solution of these 
RG-equations includes complicated summations of diagrams which are similar to solutions of a 
coupled set of Schwinger- Dyson equations. The heat kernel regularization method preserves 
the symmetries of the theory and keeps the physical interpretation of the evolution equations 
particularly simple for phenomenological applications. 

We choose as a starting action at the UV scale the NJL Lagrangian 

S.,, = J d'x[qi^dq + G {{qqf + {qifl.qf]] (1) 

which is embedded into a linear cr-model with = 1, a vanishing wave function renormalization 
constant = and an effective potential f/(<l'^) = h^^. We assume the linear cr-model to 
be a vahd description of nature below scales of 1.5 GeV. In this region the gluon degrees of 
freedom are supposed to be already frozen out. Confinement is not dealt with correctly. The 
physics is governed by chiral symmetry, which is spontaneously broken in the infra-red. The 
Euclidian action of the linear a-model in generalized form is given by 

S[^, q,q] = I d^x (^Zg{^'')qjdq + (d^^^ {d^.^') + 9{^^)mq + • (2) 

The quark and meson fields are g, q and $, where $ = (a, tt) is the 0(4)- and M = a + if'fc'j^ 
the chiral SU{2)i ® S'[/(2)R-representation of the meson fields. We consider this action to be 
preserved during the evolution and follow the generalized couplings and the effective potential 
during evolution by calculating the effective action Seff in a one- loop approximation. 

S,ff[^, q, q] = S[<^, q, q] - ^Tr log S^^ + ^Tr log {S^.^, - 23^.^3^^' S^^,) , (3) 

with 

5'3[^,q,q] 



3^i^j 

{x, y) 



6q{x)5q{y) 
S'3[^,q,q] 



6^^{x)6^i{y) 



(4) 
(5) 
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The first logarithm results from the fermion loop fluctuations, whereas the two terms in the 
second logarithm are the contributions of the bosonic and the mixed loop, respectively. The 
derivative of S^ff with respect to the evolution scale can be computed after a cut-off function 
is introduced. For this purpose we represent the fluctuation determinants by Schwinger proper 
time integrals with the cut-off function /(A;^r) = (^1 + k'^r + |(A;^r)^^ e'^^"^ 

Tr log{A) = -Tr / — e-"^/(A;V). (6) 

Jo T 

The heat kernel cut-off function / suppresses all fluctuations with momenta below the cut- 
off scale k. Going to A; ^ means to include more and more infra-red modes. Finally, all 
modes are included, because of the limiting behaviour fk^o — > 1. The ultraviolet region is left 
undisturbed due to the extra r and terms in front of the exponential. Therefore we can 
extend the r-integration to the interval [0, oo]. We further replace the masses and couplings of 
the classical action by the running masses and running couplings of the effective action Seff- 
Generically, we obtain the following evolution equation for Sejf which contains its second order 
derivatives S'^^j indicated in Eqs.(4,5): 

^--iTr/^e"<//(^)|-/(^V). (7) 



dk 2 J T dk' 

This replacement of S" by S'^jj turns the one-loop equation into a renormalization group 
improved flow equation, which goes beyond the standard one-loop renormalization group run- 
ning and includes higher loop terms successively into the proper time integral. Therefore it is 
capable to include nonperturbative physics in the strong coupling region. The evolution equa- 
tions are obtained by comparing the derivative on Seff{k) with respect to k with the formal 
expression of Eq. (2) which contains the generalized couplings and the effective potential as 
functions of k. The local potential terms are easy to evaluate and lead to nonlinear partial 
differential equations containing derivatives of the potential with respect to The evolution 
of the other couplings can be obtained in a derivative expansion which is an expansion in deriva- 
tives (momentum) over mass. This expansion is well controlled because of two reasons. Firstly, 
we approach the infra-red, so momenta become always smaller, we need only few terms in the 
expansion. Secondly, the infra-red cut-off scale k acts like a mass cut-off, so the expansion is 
well defined even if the excitations themselves are massless, like the quarks in the ultraviolet 
or the mesons in the infra-red region. 

Let me give an overview of the results we have achieved so far. We have solved in local po- 
tential approximation the vacuum evolution, the finite temperature evolution and finite baryon 
density evolution There are strong similarities between the vacuum evolution as a 



function of resolution scale and as a function of temperature. Both show the transition from 
partons to constituent quarks which has also been phenomenologically checked in deep inelastic 
scattering at fixed energy y/s as a function of photon virtuality see Ref. [0 and Fig. (2). 



Numerically, the transition temperature 27rTc is about equal to the chiral symmetry breaking 
scale k^sf, where the quarks condense in the vacuum. Going from the ultraviolet to the infra-red 
the strong attractive quark-antiquark interaction leads to a condensation of quark pairs in the 
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vacuum. At the resolution k^sb the effective mass parameter of the mesons equals zero. The 
critical exponents agree with the 0(4) critical exponents, but differ from the mean field values 



In large N^. approximation we have shown that the renormalization group flow equations are 



identical to the gap equations of the NJL model |]T5[. Further it is possible to calculate the 



eigenmode spectrum of the Euclidean Dirac operator [|T^ which previously was only known for 
very small eigenvalues in random matrix theories. 

The full equations possess extra complexity via the expansion of all running couplings and 
wavefunction renormalization parameters in As mentioned above the global form of the 
meson potential changes from the parabolic shape in the ultraviolet to a mexican hat 

shape in the infra-red. This evolution leads to a change of the minimum of the effective potential 
with resolution k. We track all couplings in Taylor series around this running minimum with 
terms up to 5th order in This leads to 20 coupled non-linear differential equations. Including 
meson loop terms the infra-red properties of the system depend only on the mass term in 
the ultraviolet and no longer on the starting ultraviolet cut-off. We have a real fixed point 
behaviour for the pion decay constant, meson quark coupling, constituent quark mass and 
quark condensate. In comparison to NJL the flow equations have predictive power since only 
one mass scale fixes all the other dimensionful parameters. 

For finite density the model has similar deficiences as the conventional NJL model. It 
overbinds and restores the chiral symmetry too early. We have improved the model by including 
cu-repulsion and plan to do this calculation including loop corrections. 

The very successful chiral perturbation theory alone cannot connect to QCD in the ultraviolet. 
Typically it starts to fail when the resolution scale is of the order of half the p mass. Hybrid 
models like the one presented here containing meson degrees of freedom and quarks extend 
to larger momentum scales, which is a definite advantage. There are, however, two problems: 
Firstly, the number of mesons in the theory should include also the vector mesons. Secondly, 
the size of the scalar Yukawa couplings increases dramatically in the ultraviolet, therefore it 
seems difficult to cross over to the asymptotically free gauge theory of QCD. For nuclei a theory 
with quarks and nucleons [T^ seems to be the most efficient way to describe the transition from 
purely nucleonic matter to quark matter at high baryon density . Our model like most of the 
other hybrid models lacks a field theoretic mechanism to avoid the appearance of free quarks as 
asymptotic states. A major progress would be to capture the smooth infra-red limit of lattice 
QCD in a field theoretic continuum picture. 
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